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>v Abstract 

a ■ 

The relativistic motion of an isolated two-body system (bound or unbound) 
of given lab energy K° in QED is separated into cms motion and relative 
motion. The relative motion equation /C L ?/> L (r L ) = contains the momentum 
j> ' eigenvalue K of the cms motion. It is greatly simplified by a binary boost to 

the atomic rest frame, where K° and K appear only in a Lorentz-invariant 
^ . combination. This boost is not a product of single-particle boosts, which are 

useful only for perturbative interactions. CPT-invariance is demonstrated, 
\ and orthogonality relations are derived. 

on 



a: 



I. INTRODUCTION 



OS. 

This paper is about the relativistic two-body problem in QED. It treats the separation of 
the cms motion from the relative motion in analogy with the nonrelativistic separation. The 
case of two spinless particles has been treated previously |] but has no practical applications. 
^ ■ When one of the two particles is a lepton (e, /i or r) and the other is spinless, the two-body 
wave function contains an additional boost: 

A = ( 7 + 75 K<x 1 /£) 1/2 , J = K°/E, K° 2 -K 2 = E 2 (1) 

in units H — c — 1. Here = (K°,~K) are the eigenvalues of P M = (idt,P) with 
P = pi + p2 ; cr\ are the Pauli matrices of particle 1, which are normally combined with the 
Dirac matrix 75 into ot\ = j^cri- The binary boost A resembles the boost Ai for a single 
free particle: 

A 1 = ( 7l + 75 K l( T 1 /m 1 ) 1/2 , ll = K° 1 /m 1 , Kf - K 2 = m\. (2) 

(An alternative form of A x which avoids <Ji under the square root is (2m 1 )" 1,/2 (K® + mi) 
(K® + mi + CKiKi), but @ is in fact more elegant.) These points will be elaborated in 
section fTJ. They are of interest for muonic helium (fi~a), mu-pionium (/i~7r + ) and also for 
precise recoil corrections in ordinary atoms. 
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Binary boosts are simplified by taking the z-axis along the total momentum K , c^K = 

(J\ Z K , and by abbreviating K/E = jv/c as K: A = Oy + ^Kai^j . The symbol (3 will 
denote the parity Dirac matrix which will be taken in diagonal form: 

^ = 70 = ( o -°i ) ' 75= (fj)- ® 

The boost which we find for two spinor particles contains a z = a lz + a 2z as well as a 
/^-dependent mass ratio 

B=(l + ^K 2 a 2 z + iipKa^ K = K/E , (4) 



fxp = M_/M + , M±=m 2 ± Pm 1 . (5) 
As a z has only the two eigenvalues 4 and 0, an alternative form of B is 

5 = 7 + npKa t <y s /2 , 7 = \fl + K 2 <r 2 j4: . (6) 

Evidently, 7 = 7 for cr 2 = 4 and 7 = 1 for cr 2 = 0. 5 contains two Pauli matrices but only 
one /3 and one 75 , i.e. B is an 8 x 8 matrix. For comparison, the two-body boost B12 used 
so far in perturbative QED is a 16 x 16-matrix, 

,812 = A X A 2 = (71 + piai/mi) 1/2 (72 + p 2 «2/m2) 1/2 • (7) 

It requires two separate conserved energies, K® and K®, momentum operators pi and p 2 
and Dirac matrices /3i, j3 2 , 751, 752- The main advantage of B is that it contains only the 
total momentum K instead of the individual momenta. It Lorentz transforms the asymp- 
totic part of the cms spinorial wave function, the radial part of which is e tkr in spherical 
coordinates (bound states have imaginary k = in; in QED, the ground state has kq = a^ 1 , 
as being the relativistic Bohr radius). The boost Bi 2 , on the other hand, requires asymp- 
totic wave functions e * k i r i e lk 2 r 2 anc [ then by analogy also two different times, e~ lK ° fl e~ lK ° t2 . 
The derivation of B is given in the next section. As in the nonrelativistic case, the shape 
of the potential V(r L ) is irrelevant in this derivation (r L — r\ — r 2 ). Anomalous magnetic 
moments can be included in the cms equation M but are as yet excluded from the boost. 
Thus the present treatment applies to the bound states of two charged leptons, which may 
be called leptonium (muonium e~fi + , antimuonium e + /i _ , positronium e~e + ), as well as to 
their scattering states up to K° — > 00 and / or E — ► 00. 

The relativistic separation of the cms coordinate R reads |l| 

r 1 = R+^(l- Am 2 /£ 2 )r L , r 2 = R - ]- (l + Am 2 /E 2 ) r L (8) 

with Am 2 = m\ — m 2 = —M + M_. For loosely bound states, E 2 w (m\ + m 2 ) 2 reduces (||) 
to the nonrelativistic transformation except for the index L in r L which still signifies the lab 
system. This additional index is necessary because r L is Lorentz-contracted in the direction 
K. The uncontracted cms variable will be denoted by r*. In the following, the z-axis is 
taken along K, such that x L and y L are not Lorentz-contracted, and 



2 



£ L = z*/j, p Lz = 7p*. (9) 

The usual cms potential V* = —a/r* (a = e 2 ) which is abstracted from the QED Born 
approximation is Lorentz-contracted in the lab system. In the 16-component Dirac-Breit 
equation which is applied successfully in atomic theory, this contraction is not evident, but 
a careful analysis shows that it is represented by the Breit operators [0. With a simple 
V(r) = —a/r, the cms equation is free of Breit operators. However, this form in fact 
requires a third variable transformation which affects only the distance, r — > r*. In other 
words, the Fourier transform of the first Born approximation of the cms scattering amplitude 
produces a potential V(r) where r is a quasidistance from the Dirac-Breit point of view. 
These procedures are up to now only approximately Lorentz-invariant and are not repeated 
here. We merely wish to reserve the index-free r for the quasiposition, which may or may 
not coincide with r*. 

An important application of this 8-component formalism is to positronium. The relation 
m,2 = mi simplifies (||) to ri = R + r L /2, r 2 = R — r L /2. In this case, B can be used in the 
small components (/3 = — 1 , /ig = oo) only for a 2 z = , in which case the large components 
have a 2 = 4 . The opposite case requires a 75-transformation, which will be explained in 
section [D]. 

We have found in the literature a formula which resembles (f|): Replacing E — > m\ and 
Hp — > 1 , (f|) becomes identical with a boost for a single free particle of spin 1 || . 

II. DERIVATION OF THE 8 COMPONENT EQUATION 

A convenient starting point in the lab system is at present the 16-component Dirac-Breit 
equation 

7r£ - mift - maft - PbWdl = , tt° = id t - V 12 (10) 



where P B is the Dirac-Breit momentum operator |l[], and the indices D and L stand for 
Dirac and lab system, respectively. It contains pi and P2 multiplied by Breit corrections 
1 — V12 {ctiCx-2 + Ci\ r oi2r) /2, which contain low-energy approximations. In the future, one 
would like to formulate the Feynman rules for the fermion-fermion scattering amplitudes 
in terms of 8-component spinors, both in the lab system and in the cms, and then derive 
the interactions via Fourier transformations. At present, this derivation is incomplete even 
in the cms. The cms equation has been abstracted from (fLOf) , but most of its scattering 
amplitudes have been checked against the Born approximation for arbitrary energies. For 
the present construction of the free two-particle boost, we only need three properties of 
(H0) : (i) ft i s translationally invariant in space and time, 

id t 1pT3L = A'VdL , (Pi + P2) ^DL = K^ DL , (11) 

(ii) it is exact at all energies in the asymptotic region V12 = and (iii) it contains only a 
single time t, which contributes a factor e~ tK 1 for stationary states. The resulting time- 
independent Dirac-Breit equation is again (10), but with 7r° = K° — V\2- It contains the six 



variables Ti and r 2 , precisely as in the nonrelativistic case. For V\ 2 = 0, we thus have 

\K° - rmPi - m 2 (3 2 - Pie*i - P2CK2) V^dl = . (12) 
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In the parity basis in which (3\ and fa are diagonal, the components of ip DL may be labelled 
by the index g for a large component and / for a small component of either particle: ip gg , ip g f, 
ipfg and ipff. In the following, these components are rearranged into sums and differences 
as follows: 

-(^\- 1 (^SS+^fA v . -(Xg\_ 1 (ipgg-'&fA /ion 

The index D stands for Dirac, the index L for the lab system. The factor 2 -1 / 2 makes the 
transformation unitary The 16-component equation thus assumes the form of two coupled 
8-component equations, 

(K° - 75 p+) ^ DL = (3M + Xdl , (K° + 75 p_) Xdl = ^ DL , (14) 



P± = PiO"i ± p 2 cr 2 , (15) 

with /3,75 and M + defined in ([3]) and (|). In words, 75 exchanges the single index g with 
the single index / in ip and %■> while f3 multiplies / by — 1. The original matrices Pi, P2, 751 
and 752 are now obsolete. Elimination of Xdl by means of the first equation (|T4]) yields 

(if + 75 p_) /3M7J 1 (V - 75 p + ) tp DL = (3M + tfj DL . (16) 
Upon multiplication by (3M +) this becomes 

(V? _ M 2 _ KOp^-i^ _ ^o p+7s + = 0. (17) 

The denominator M_ = m 2 — /3toi appears here because of 75/? = —P75. At this point, we 
replace ri and r 2 by R and r L according to (§), which leads to 

Pl = p L + i(l + Am 2 /£ 2 )K, p 2 = -p L + ~ (l - Am 2 /E 2 ) K (18) 

where ([□]) has already been used, and E 2 = K° 2 — K 2 as in ([I]). ip Dh is now a function of 
the vector r L and of the parameters K° and K (with the z-axis along K). We now show 
that ip DL can be reduced to a function of the components r Lt = (x L , y L ) and z* = 7z L which 
contains K° and K only in the combination K° 2 — K 2 = E 2 . This justifies the choice of the 
transformation (§P, a posteriori. 
The combinations p± of (|Uj) become 

/ , Am 2 \ 1 

P- = PlCT* + 7P? H — K cr z H A Acr z , . , 

^ 2E 2 ) 2 ' (19) 

/ Am 2 \ 1 
p+ = p L A<x t + I 7 p* + 2^-^ J Aa, + -Ka z , 



cr = cr 1 + cr 2 , Acr = a 1 - cr 2 , cr t = (a x ,a y ). (20) 

Here we have denoted pu^t = PlC* in order to save one index. The operators pu and p* refer 
to the cms, but we avoid the notation p = —iV which is reserved for the quasimomentum. 
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In the Pauli spinor space xi X2 , cr is symmetric in 1 and 2. It transforms spin triplets into 
triplets and annihilates spin singlets. Act is antisymmetric and exchanges singlets (ip LS ) with 
triplets (ipu)- It turns out that the following separation of mass factors is useful: 



V-'e 



(21) 



The index D disappears here; the new components ipu and ip hs form an 8-component spinor 
ip h , which satisfies the following equation: 



K ° 2 - M 2 + + ii/p-p+ - K° 



2p L cr lt 4 ( 2- ik 
K 



L\m 2 K" 



E 2 



+ — (V/Ao^ + iipa^ 75}^ = 0. 



Here we have used cr + Act = 2a \. The expression p-.p + simplifies considerably, 

Am 2 



P-P^ 



E 2 



-K z + 2K 1P * 



(22) 



(23) 



Remembering Am 2 = -M . M . the form M + Am 2 K 2 /E 2 M- = —M\K 2 jE 2 is combined 
with — M\ into —M\(\ + K 2 / E 2 ) = — M+7 2 , 7 = K°/E. One may thus remove one factor 
7 from (P2"p. The result is JC L ip L = 0, where 



(E - Mi ) 7 + 2KnpY z - 2Ep L cr lt + (2Ejp* z - M + M_K)a lz + 



+-KE 2 (np 1 Acr z + npcr t 



75 



(24) 



with K = K/E as before. In the cms, i.e. K — 0, 7 = 1, this reduces to 



£V = , K = E 2 - Ml - 2E(p L <r lt + p> lz b 
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The operator /C is independent of cr 2 and well-behaved for mi = m 2 , where M 2 
annihilates the small components. It remains to find a transformation ip h 
reduces /C L to /C after multiplication by another matrix B from the left 



(25) 

2m 2 {\+(3) 
Bip which 



(26) 



BJC h B^ = JCtp = 0, K = B)C L B . 
As the operator — 2£'p L <xi t occurs in (|25|) in the same form as in (24[), one needs 

B(T u = <r lt B-\ (27) 
The form (|2~4f) suggests that may not contain er lt and cr 2 f. Next, one may note that 

<y\ z cr\t = —cr\tO\ z , Aa z a- U = —cr\ t a z . (28) 



After a few fruitless attempts, one solves (^) with 
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B = (l + -K 2 a 2 + nKam) , B' 1 = (l + -K 2 a 2 - fipKa^j , (29) 

B=(l + \k 2 A<j 2 z + ^KAa zll ^j V2 . (30) 

It remains to show how the eigenvalue oo of /ig (|J) is avoided for positronium. When the 
large components have a 2 = 4 , a z = ±2 implies mi = rrij ± 1 (mi and mj are the eigenvalues 
of L z and J z = L z + a z /2). These states have parity (— (J 2 = j(j + 1)). The orbital 
parity of the small components is opposite, i.e. (— The spin function can be either singlet 
or triplet, and in either case mi = mj ensures a z = , i.e. the combination {ipo z can be taken 
to vanish in (^). Thus in the application to positronium, one must discuss the spin structure 
for mi 7^ m 2 and then take mi = m 2 in the final forms. When the large components have 
a z = 0, a chiral transformation ip DL = 75^ DLjC / l transforms fip into /il , such that the above 
argument applies again. 

III. THE KLEIN DIRAC BOOST 

The 4-component Dirac-Klein-Gordon-Breit equation for systems such as yU~7r + reads 

1 

(Vf + Kl _ £ 2 _ 2m 1 7r^) < = 75 ({Pi, ?r L Vi + 6i) ^ (31) 
where &i is the Breit modification of the momentum operator, 

h = -V U (<TlP2 + OXrVlr) , (32) 

and JCi = m\ +p 2 . The coordinate transformations @ and @ give 

/Ci - /C 2 = 7 2 Am 2 + 27p*X . (33) 
Setting now Vi 2 = and extracting one factor 7 , one obtains /QVl = > with 

= (E 2 + Am 2 ) (7 - j 5 Ka lz ) + 2Kp* z - 2 l5 E ( 1P * z a u + p L< r lt ) - 2m x (3E . (34) 
Writing the lab spinor %p h as a boost A times the cms spinor ip' , ip'^ = Aip' , one finds 

/Cty' = 0, K.' = AJC' h A, A=( 1 + l5 Ka lz ) 1/2 , (35) 

/C' = E 2 + Am 2 - 2£mi/3 - 2E (p* z a z + p L a t ) . (36) 

The asymptotic energy E\ of the spinor particle is E\ = (E 2 + Am 2 )/2E , such that the first 
two terms in (|36|) arise from 2EE\. As /C' has the Dirac operator structure, it is clear that 
its eigenvalues depend on E 2 + Am 2 and —2Emif3 only via (E 2 + Am 2 ) 2 — Am\E 2 . Setting 
now 

V/ = (2k 2 Ey 1/2 (E x - (3m,) 1 ' 2 (E 2 - M 2 f /2 ^ , k 2 = E 2 — m\ (37) 
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one finds that JC' is transformed into JC as given in ( p5|) (the resulting boost is rather 
different, however). The point k 2 = comprises two thresholds at E 2 = (m% +m 2 ) 2 and two 
at E 2 = (mi — fn 2 ) 2 . In our example of a /i~7r + system, E = m\ + m 2 is the fi~ir + threshold, 
E = —rrii — m 2 the fi + ir~ threshold, E = m 2 — nil the fi + n + threshold and E = m% — m 2 
the yU _ 7r _ threshold. For mi = m 2 , one factor E can be separated from fl36|), leading to 
JC' — E — 2mi(3 — 2{p*a z + p L a t ) . This completely eliminates the doubly-charged channels, 
which now require a separate equation. A similar decoupling occurs in the double-Dirac 
case for mi = m 2 , which is discussed below. 



IV. ORTHOGONALITY RELATIONS 

A single 4-component Dirac spinor has components ip R and ip L in the chiral basis, where 
7 5 is diagonal. f3 and 75 simply change their places in (f|), i.e. (3 exchanges ip R and ip L . As 
(3 is part of the parity transformation, ip R and ip h are not parity eigenstates. They do form 
separate representations of the Lorentz group, however. 

In the present case, both ip and x can have /3 diagonal and remain separate under Lorentz 
transformations. Elimination of ^ DL from (|T4|) leads to an equation for x DL in which the two 
brackets of fljjD are interchanged, which is equivalent to the substitution p + < — > — p_ . The 
equation corresponding to (|ITD for Xdl is thus 



( K ° 2 - M 2 + + (flf P+ + p_) K° l5 + fx/p+p-) Xdl = . (38) 
It will be shown in the following that the orthogonality relations require both ip and \. The 



mass separation analogous to (]2l|) is 



^U&J- (39) 

The relation analogous to fl2"4"|) is /C XL x L = , with 

/C XL = (£ 2 - M 2 ) 7 + 2^^: - [2£p L o-ii + (2^ - M+M-K)<r lz + 

+ \KE 2 (npAa z + ^V 2 ) ] 75 , (40) 

The boost is now different, x L = B x x, where 

B x =(l + \k 2 A<j 2 z + ^tfAoyns) ^ = £ f (41) 

where B has been given in (^). The hermitian conjugation of B in (^) merely exchanges the 

positions of fi^ 1 and 75. Observing (/XgVs) = 75/U3 1 = yU/375 , one readily verifies (|4]). The 
coordinate transformation (^j) implies d 3 rx d 3 r 2 = d 3 r L ti 3 i? , and the orthogonality of states 
with different momenta K follows simply from / d 3 R exp{— z(K — K')R} = (2tt) 3 5(K — K'). 
A subsequent boost allows one to discuss the orthogonality of the functions ip(r L ) in the 
cms. The relevant equations here are JCip = and JC\ — , with JC given by (|25f) . The latter 
equation follows from JC XL x L = ior K = , j = 1 , i.e. JC X = JC in the cms. 
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In the ordinary Dirac equation, neither the potential V(r) nor the momentum p = — «V 
appear in the orthogonality relations. Both properties are also achieved here, but the details 
are surprising. The general operator form of /C will be needed for all values of r, not just for 
r — > oo where V(r) = 0. This form can be derived from the Fourier transform of the QED 
Born approximation for elastic scattering 1 + 2 — > 1' + 2' in the cms: 

K = E 2 - Ml - 2Ecxp - 2EV(r) - AV'(r) , (42) 

ol = 7 5 ct! , A = %a r a 1 rr 2 = (a x cr 2 ) r + i^ 5 a 2r , (43) 

and V = dVj dr , a r = err as usual. The derivation from the Dirac-Breit equation leads 
to a more complicated form in the variable r* = (z*,x h ,y h ) , which involves Breit operators 
and also V 2 (r*). However, it turns out that the substitution 

r*^r + a/2E (44) 

does reduce the Dirac-Breit expression to ([O) at low energies. As the QED Born approxima- 
tion is free of low-energy approximations, ([42]) is much more practical than the Dirac-Breit 
equation. 

To derive the orthogonality relations for the solution of a Hamiltonian equation Hip = 
Eip , one writes Hipj = Ejipj , {Hipi)^ = Eiip\ , multiplies the first equation by tpj , the 
second one by ipj , subtracts the second product from the first one, and integrates over all 
configuration space: (E t - Ej) J ipj^j = 0. The method can also be applied to the Klein- 
Gordon equation, \{E — V(r) ) 2 + V 2 — m 2 ]ip = (in units Ti = c = 1), but the square of 
E — V produces a weight Wij = Ei + Ej — 2V, i.e. / ip* Wij ipj = 5ij . The case at hand is more 
complicated because A is not hermitian. Its hermitian part is the recoil-corrected hyperfine 
operator. Its antihermitian part has zero expectation values in fine structure eigenstates; 
it is needed in positronium where fine and hyperfine structures are comparable (the energy 
eigenvalues remain real). Fortunately, the equation K, x x = has A replaced by A' : 

}C X = E 2 - Ml - 2Ectp - 2EV(r) - AV'(r) . (45) 

One may thus envisage orthogonality relations / xl ^ij i>j — $ij , where the antihermitian 
component disappears with A^ — A = 0. However, the factor E multiplies not only V(r) as 
in the Klein-Gordon equation, but also ap. To avoid p in the orthogonality relations, E 
must be divided off. But then A/E is both non-hermitian and energy-dependent, in which 
case V will remain in the orthogonality relations. 

For V = —a/r, one may introduce a dimensionless scaled variable, 

r = Er d/dr = E^d/dr p = p/E, (46) 

and divide /C by E 2 = s: 

(2ap + 2V(r) + AV'(r) - 1 + M 2 + /s) ip(r) = 0. (47) 

Using the corresponding equation for x ■> one arrives at 
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1 - sT 1 ) / X \ Ml ^ d 3 r = , s t = El , (48) 
xlM 2 + i] 3 d 3 r = 6 lJ . (49) 



Remembering M\ = m\ + m\ + 2mi?7i2/3 , this is a simple generalization of the static limit 
m\jmi = 0. For positronium, the small components do not contribute to fl49|). 

Equation (|48| ) is explicitly CPT-invariant: Every bound state Sj has two different eigen- 
values Ei , namely Ei = ^/Si = m Ai and Ei = —■sfsi = —m Ai , where m Ai denotes the atomic 
mass in the state i (an excited atom is heavier than its ground state). Returning now to 
the time-dependent cms equation with idt = E , one finds that Ei = —m A i belongs to the 
antiatom of mass m A i , i.e. m A i = m A i . Positronium is its own antiatom, of course. 

Vacuum polarization introduces an extra scale into V(r). The equation remains CPT- 
invariant, but (^) becomes more complicated. In the Klein-Dirac case of section [HI] , the 
spinless particle normally has an extended charge distribution, leading to V(r) ^ —a/r for 
small r. The simplest way out is of course to take V(r) as an arbitrary function of f in 
(ff7|). In any case, finite charge distribution models need some tuning in order to conform 



with CPT. 

It should be stated that the range of the dimensionless radial variable f is < f < oo 
both for atoms and for antiatoms. In the old variable r, antiatoms have negative distances. 
Of course, this makes as little sense as the claim that antiatoms By backwards in time. 

We conclude with two comments on possible applications of the new equations: A vector 
potential A(r, t) is included in (|T^) by replacing pj — > iZi = pj+gjA(rj, t). This allows one to 
calculate relativistic recoil effects in positronium de-excitation, or in the e + e~-recombination 
into positronium. In addition, the Lamb shift may be calculated from these processes via a 
dispersion integral. 

The second comment concerns the case mi = rri2 . Here M\j s vanishes for the small 
components ipj of ip . Calling for brevity 

V = V(r), n = p+ l -VV<T 1 <r 2 , (50) 

one has ipf = (1 — 2V)~ 1 2(Ti7rip g . The large components satisfy the equation 

~4mj/s -1 + 2V + 4cr 1 7r (1 - 2V)~ l (T i7T J ip g = . (51) 

For V = + a/r , it has only scattering states, which would be e~e~ and e + e + in the case of 
two electrons. If it were possible to extend the present formalism to an external Coulomb 
potential, one would arrive at a theory which isolates the e~e~- and e + e + - channels from 
the e~e + -channel. This would dispense with positive-energy projectors and simplify the 
relativistic variational calculation of atomic ground states. 
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